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Abstract 


An algebraic one-dimensional model on the membrane-electrode-assembly (MEA) of direct methanol fuel cell (DMFC) is proposed. Non-linear 
regression procedure was imposed on the model to retrieve important parameters: solid polymer electrolyte conductivity Km, exchange current 
density of methanol electro-oxidation at anode catalyst surface igq rer, and mass diffusivity of methanol in aqueous phase within the porous electrode 
D, that correspond to the experimentally measured polarization curves. Although numerical iteration is required for a complete solution, the explicit 
relationships of methanol concentration, methanol crossover rate, oxygen concentration and cell discharge current density do provide a clear picture 
of the mass transport and electrochemical kinetics within the various porous media in the MEA. It is shown the cathode mixed potential induced 
by the parallel reactions of oxygen reduction and oxidation of crossover methanol elucidates the potential drop of the cathode and the decrease 
of the cell open circuit voltage (OCV). Methanol transport in the membrane is described by the diffusion, electro-osmosis, and pressure induced 
convection. Detailed accounts of the effects of anode methanol and cathode oxygen feed concentrations on the cell discharge performance are 


given with correlation to the physical structure and chemical compositions of the catalyst layers (CLs). 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Direct methanol fuel cell is the recent focus of power sources 
on portable electronic devices such as laptop computer, cellu- 
lar phone, and video/audio player, to name a few [1]. Although 
direct methanol fuel cell (DMFC) has the benefit of small system 
size and weight, which can be operated at low temperature, it 
is suffered from the slow electrochemical kinetics of methanol 
oxidation at the anode and the voltage loss due to methanol 
crossover through the membrane separator [2—4]. In addition, 
water management is another crucial issue for a compact design 
of portable DMFC [5-8]. Wang and co-workers [7,8] indicated 
that a severe water loss occurs in the anode due to the electro- 
osmotic drag through the polymer electrolyte membrane, which 
prohibits the use of concentrated methanol fuel so that addi- 
tional equipments are required to replenish water in the anode 
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from that in the cathode. Further, the conceived slow oxygen 
reduction in the cathode could worsen the cell performance 
in many occasions [9]. Experimental efforts have been carried 
out to search and formulate more efficient, carbon monoxide 
resistant anode catalysts [10,11], and to modify the solid poly- 
mer electrolyte membrane for a better methanol/water barrier 
[12,13]. All the materials relevant researches were devoted to 
the enhancement of single cell performance, and with an eye 
on a favorable chemical/physical structure of the membrane- 
electrode-assembly (MEA) [14,15]. 

Many mathematical models were developed to assist the 
understanding and improvement of MEA fabrication technique 
on DMFC [16-29]. A comprehensive account on the current 
advancement of modeling methods and the associated exper- 
imental diagnostic techniques on the polymer electrolyte fuel 
cells was given by Wang’s review article [16], in which a 
large body of relevant literature was documented. Scott and 
co-workers [17—21] focused on the mass transport of reac- 
tants through the porous electrodes and the interplay with 
associated electrochemical reactions. Their simulation analyzed 
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Nomenclature 


da 


ac 


i oM, ref 


1o0,ref 


active area per unit volume in anode catalyst layer 
(cm!) 

active area per unit volume in cathode catalyst 
layer (cm7!) 

active area per unit mass of PtRu catalyst 
(cm? g7!) 

active area per unit mass of Pt catalyst (cm? g7!) 
anode channel concentration of methanol 
(mol cm~?) 

proton concentration 
(mol cm~) 

methanol concentration (mol cm7?) 

reference methanol concentration (mol cm~?) 
oxygen concentration in cathode catalyst layer 
(mol cm~?) 

oxygen gas concentration at cathode catalyst 
layer/gas diffuser (mol cm7?) 
cathode channel oxygen gas 
(mol cm~?) 

reference oxygen concentration (mol cm~?) 
methanol concentration in anode catalyst layer at 
position z; (mol cm~?) 

methanol concentration in anode backing layer at 
position z; (mol cm~?) 

methanol concentration in anode catalyst layer at 
position zy (mol cm~?) 

methanol concentration in membrane at position 
zi (mol cm~?) 

methanol concentration in membrane at position 
zn (mol cm™°) 

methanol diffusivity in anode catalyst layer 
(cm? s7!) 

effective diffusivity of methanol in anode catalyst 
layer (cm? s7!) 

methanol diffusivity in anode backing layer 
(cm? s7!) 

effective diffusivity of methanol in anode backing 
layer (cm? s7!) 

oxygen diffusivity in cathode catalyst layer 
(cm? s7!) 

effective oxygen diffusivity in cathode catalyst 
layer (cm? s7!) 

oxygen gas diffusivity in cathode gas diffuser 
(cm? s7!) 

effective oxygen gas diffusivity in cathode gas 
diffuser (cm? s7!) 

proton diffusivity in membrane (cm? s7!) 
methanol diffusivity in membrane (cm? s~!) 
Faraday constant (96,500 C mol7!) 

reference methanol oxidation exchange current 
density (A cm~?) 

reference oxygen reduction exchange current 
density (A cm~?) 


in membrane phase 


concentration 


Po), ref 
Po), ref 


Veell 


cell discharge current density (A cm7?) 
cross-over current density (A cm~?) 

oxygen reduction current density in cathode cat- 
alyst layer (A cm~?) 

local transfer current density within catalyst layer 
(A cm~?) 

potential dependent rate constant of methanol oxi- 
dation (s7!) 

potential dependent rate constant of oxygen 
reduction (s7!) 

permeation constant of pressure induced convec- 
tion (cm? s~! atm!) 

solubility constant of oxygen in liquid phase 
equilibrium constant of methanol between anode 
backing/anode catalyst layer 

equilibrium constant of methanol between anode 
catalyst layer/membrane 

anode catalyst layer thickness (cm) 

anode backing layer thickness (cm) 

cathode catalyst layer thickness (cm) 

cathode gas diffuser thickness (cm) 

membrane thickness (cm) 

cathode Pt catalyst loading (g cm~?) 

anode PtRu catalyst loading (g cm7?) 

molecular weight of O2 (gmol~!) 

molecular weight of N2 (gmol7!) 

molecular weight of H20 (g mol~!) 

electrons transferred of methanol oxidation (=6) 
electrons transferred of oxygen reduction (=4) 
methanol flux (mol cm~? s7!) 

anode pressure (atm) 

cathode pressure (atm) 

critical pressure of O2 (atm) 

critical pressure of N? (atm) 

reference pressure of CO2 (=1 atm) 

reference pressure of O2 (=1 atm) 

universal gas constant (8.314J mol~! K~!) 
universal gas constant 
(82.06 atm cm? mol ~! K7!) 

absolute temperature (K) 

critical temperature of O2 (K) 

critical temperature of N2 (K) 

reference methanol oxidation open circuit voltage 
(V) 

standard potential of methanol oxidation (V) 
reference oxygen reduction open circuit voltage 
(V) 

standard potential of oxygen reduction (V) 
anode electrode potential (V) 

molar volume of O2 at normal boiling point 
(=14.8 cm? mol!) 

cathode electrode potential (V) 

cell voltage (V) 
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Greek letters 

a anodic transfer coefficient of methanol oxidation 
Oc cathodic transfer coefficient of oxygen reduction 
Ea void fraction in anode catalyst layer 

Eb void fraction in anode backing layer 

Ec void fraction in cathode catalyst layer 

Ed void fraction in cathode gas diffuser 

Da polymer electrolyte potential in anode (V) 

Øg association parameter of H20 (=2.26) 

De polymer electrolyte potential in cathode (V) 

Na electrode over-potential in anode (V) 

Nc electrode over-potential in cathode (V) 

Km proton conductivity in membrane phase (S cm7!) 
LB viscosity of water (=1.45 cP) 

pc density of carbon (g cm~?) 

pPt density of Pt (g cm~’) 

PRu density of Ru (g cm7?) 

a electro-osmotic drag coefficient (cm? mol~!) 


individual effects of anode and cathode compartments, and the 
diffusion coefficient of methanol in the electrode matrix was 
retrieved in correlation with the experimental data [17,19,20]. 
Although methanol crossover phenomenon was stressed in their 
model, an argument of “poisoning effect” on cathode was 
employed to explain for the cathode voltage drop caused by 
crossover. That is, active area within the cathode catalyst layer 
(CL) is blocked by the adsorbed methanol that crossed, and 
polarization of oxygen reduction is increased due to the less 
available active area. However, such account was not supported 
by the physical observations. Wang et al. [30] measured the 
chemical compositions of cathode effluent of a DMFC and 
indicated that the crossover methanol is oxidized to CO2 com- 
pletely. Ren et al. [31] found the methanol permeation rate is 
inversely proportional to the membrane thickness at a given 
cell current density, suggesting the crossover is actually dif- 
fusion limited. These experimental findings showed that the 
over-potential imposed on the methanol oxidation at cathode 
is large so that methanol is oxidized immediately when reach- 
ing the cathode catalyst layer. Further, the intermediate poison 
species carbon monoxide will not exist at such high over- 
voltages at the cathode as indicated by the electrocatalytic 
studies of Dinh et al. [32], which is justified by the theo- 
retical prediction of cathode voltage response in the present 
work. 

Mixed potential concept was adopted on several theoretical 
modeling of DMFC [4,22—24,26—29]. The methanol crossover 
causes an internal oxidation current in the cathode catalyst layer, 
reduces the available cathode potential, and drags a decline of the 
cell voltage. Elaborate finite difference formulation were used 
by Zhang and Wang [24] and Guo and Ma [26], which is not 
computationally efficient and also difficult to be implemented 
for the in situ data analysis of actual fuel cell system. Wang and 
Wang [27] did consider the complex two-phase flow in channels 
using computational fluid dynamics technique; however, both 


the anode and cathode catalyst layers were simplified as planes 
of no thickness. Such simplification was later improved by the 
3D model of Liu and Wang [28,29], in which the effect of non- 
linear electrochemical kinetics within the finite catalyst layers 
was included and detailed water crossover effect was elucidated. 
The simplified algebraic formulation by Garcia et al. [23] has a 
neat closed form for the methanol concentration distribution but 
a linear profile had to be assumed in the anode catalyst layer in 
their derivation. In view of the non-linear behavior of electro- 
chemical kinetics, such assumption is apparently invalid at the 
condition of large discharge currents. Guo and Ma’s analytical 
approach [26] included the channel concentration consumption 
effect; however, a constant concentration was adopted for either 
methanol or oxygen within the catalyst layer. The purpose of 
the present study is to device a semi-analytical method that cor- 
rectly treats the mass transport phenomena occurring in various 
media of the MEA, as well as the relevant electrochemical kinet- 
ics within the active layers. The model includes the methanol 
crossover effect, which is efficient for the estimation of pertinent 
MEA kinetic parameters. 


2. Mathematical model 


A schematic of the membrane-electrode-assembly of DMFC 
is depicted in Fig. 1. The cathode and anode electrodes are com- 
posed of carbon supported catalyst paste that normally brushed 
onto the respective gas diffuser and backing layers. A complete 
sandwiched MEA is fabricated by hot-pressing a SPE membrane 
(for example, Nafion® membrane) in-between the cathode and 
anode electrodes. Cathode catalyst paste is composed of car- 
bon particle supported platinum catalysts, while that of anode 
catalyst paste is similarly prepared but with carbon supported 
platinum/ruthenium alloys. The introduction of ruthenium is to 


d 
cathode gas diffuser 


anode backing layer membrane 


anode catalyst layer cathode catalyst layer 


Fig. 1. Schematic of the model region of the membrane-electrode-assembly, not 
to scale. 
Reaction in anode: 


CH30H + H20 —> CO. +6Ht +6e7 


Reactions in cathode: 


O2 + 4Ht + 4e7 > 2H20 


CH30H + H20 —> CO, +6H* + 6e” 


K.-M. Yin / Journal of Power Sources 167 (2007) 420-429 423 


promote de-sorption of the intermediate CO molecule on Pt for 
the recovery of a free surface for the oxidation of methanol. 
Both anode and cathode pastes need to be impregnated with 
recast SPE for an efficient proton transport within the catalyst 
layers. Cathode or anode diffusion layer is made of either woven 
carbon cloth or carbon paper for an even distribution of oxygen 
gas or aqueous methanol. 

The modeling region does not include the fluid flow in the 
channels within the bi-polar plates. Thus, there is no concen- 
tration depletion considered in the conduits; only mass transfer 
resistance in the diffusion layer and the catalyst layer is included. 
Other assumptions are listed: 


(1) The system is isothermal and operated in the steady state 
condition. 

(2) The possible micro-convection and the increased gas void 
fraction in the anode backing layer, caused by the CO2 
generation in the anode catalyst layer, are not considered 
explicitly in the model. Nonetheless, the gas bubble effect 
will be reflected on the estimated effective mass diffusivity 
of methanol in the porous medium when fitting the model 
to actual experimental polarization data. 

(3) Water is fully flooded in the anode backing layer, anode 
catalyst layer, membrane, and the cathode catalyst layer, 
but not in the cathode gas diffuser. 

(4) The methanol crossover through membrane is oxidized 
instantaneously at the membrane/cathode catalyst layer 
interface due to the large cathode over-potential imposed 
on the methanol [23,26]. This is justified in the model pre- 
dictions. 

(5) The electronic resistance in carbon matrix is neglected. It 
should be mentioned that various internal resistances arise 
due to the different fabrication procedures adopted for the 
MEA. In addition, extra contact resistance appears when one 
assembles a single cell. Thus, the membrane conductivity 


The boundary conditions at z=0 and z= Lp are as follows: 


Cm =C atz=0 (3) 
and 
Cm = Œ@ = K\C? atz= Lp (4) 


It is assumed fast equilibrium preserved at the interface of anode 
backing layer/anode CL described by Eq. (4) [23,25]. The con- 
centration is 


KıC? — Cp 
= —— 4 


C 
M A 


+ Cb (5) 


2.2. Methanol conservation in the Nafion® membrane 


At steady state, methanol flux in the membrane is constant, 
which is the resultant of three transport mechanisms: diffusion, 
electro-osmosis, and pressure induced convection. 

dCm Teel AP 


Nu = — Dm —— "C kpC 6 
M me z +s M-p EMT (6) 


The electro-osmotic drag force is proportional to cell current 
density and local concentration of methanol as suggested by Ren 
et al. [31]. Darcy type fluid flow is dependent on the pressure 
difference between cathode and anode as well as the methanol 
concentration in the membrane [20]. Wang and co-workers [7] 
analyzed the effect of cathode hydrophobic carbon micro-layer 
on the hydraulic backing pressure. An extra capillary force pro- 
hibited the water crossover and the net water transport coefficient 
is reduced. Our model can easily be modified to include such 
feature once the pertinent thermodynamic properties are avail- 
able for this additional micro-layer. With a linear coordinate 
transformation = z — (Lp + La), Cm can be solved with the 
prescribed conditions: 


— cm 3— 
estimated in the present study should be regarded as the Sin S o 
total effects from different current conveyers of polymer and 
electrolytes and solid matrix media. 

Cm = Ch at? = Lm (8) 
This yields 
ips (Cyt — CID expl- ((kp A P/LmDm) — (5'Icen/ Dm F))2] + Cir — Cy expl—((kp AP/LmDm) — (¢/ Icett/ Dm F)) Lm] (9) 
Le exp[—((kp AP/LmDm) — (S"Tcet/ Pm F)) Lm] 
(6) Water crossover and its effect on the variation of methanol The methanol crossover rate is then given 
concentration are not considered in the present study. = (Cit — Cir )((kp AP/Lm) — (s Icen'/F)) 
~~ exp[—((kp AP/Dm) — (ç' Ie Lm/DmF))] — 1 
2.1. Mass balance in the anode backing layer Ber i P Pe 
kp AP I 
cr ( p ç =) (10) 
Mass conservation of methanol is Lm F 
Sm sü (1) Eq. (10) is simplified to 
: Ney = KICA A P/Lm) — (S Teon/ F) wo 
The diffusive flux of methanol is A A exp((kp A P/Dm) — (5'lenLm/DmF)) 
Nu = - D“ dCm (2) assuming an instantaneous equilibrium between the anode cat- 
b dz alyst layer/membrane interface Cy = KnCi [23,25], and an 
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instant consumption of methanol at the membrane/cathode cata- 
lyst layer interface Ciy ~ 0 [23,26]. The leaking current density 
caused by the oxidation of crossover methanol is simply 


Keak = nM FNM (12) 


in which ny is the number of electron transferred. 


2.3. Kinetic expression in the anode catalyst layer 


A first order electrochemical reaction on methanol is assumed 
within the anode catalyst layer: 
J 
NM F 


= kCm (13) 


where j is the local transfer current density within the anode 
active layer. The potential dependent rate constant k is defined: 


k : 1 Qana F (14) 
= al — ex 
a OMS iM PCM ret p RT 


The active area per unit volume is a4; the exchange current den- 
sity ioM.ref Corresponds to a specified reference concentration 
CM.ref- @a and na are the anodic transfer coefficient and the anode 
over-potential, respectively. na is defined as 


2.4. Mass conservation in the anode catalyst layer 
Mass conservation is derived for the methanol diffusion with 


electrochemical reaction: 


2 
eff d CM 
a dz2 


— kCm = 0 (18) 


Letting z’ = z — Lp, Eq. (18) is solved 


Cm = C? cosh a z’ | + |cosech ae, La | CA 
=U a | SI 
Dpef Deff 
k a : k 1 
— coth pet“ CI sinh Delt (19) 


under the given boundaries: 


CM = Ci atz = 0 (20) 
and 
CM = Ci atz’ = La (21) 


Since the diffusive flux of methanol at the interface of anode 
backing layer/anode active layer should be continuous, the flux 
continuity requirement at z’ = 0 is used to derive 


a (DETCy/Ly) + yk DS cosech(/k/ DSF La)Ci 
i (DË K1/ Lp) + y/kD coth(/k/ D La) 


(22) 


Na = Va — a — UM (15) Another condition required to determine Cj; is the flux continuity 
constraint at the anode catalyst layer/membrane interface (z’ = 
L,). The final result is 
4 KDD Cy / Lp) cosech(\/k/ DF La) /((De" K1/ Lb) + \/kD& coth(y/k/ DF La)) 


c= 
H (Ku((kp AP/Lm) — (¢' Teen /F)) expl—((kp AP/Dm) — (6 Leet Lm/ Dm F))])/(1 — exp[—(kp A P/ Dm) 
—(6/ Feen Lm/ Dm F)I) + (/kDS"( DE" K1/ Lo) coth(/k/ DF La) + kD / (CDR K1/ Lo) + /kDE* coth(y/k/ DF La) 


Va is the anode electrode potential, ®, the SPE potential in 
anode CL, and Uy the open circuit potential (OCP) of methanol 
oxidation at the reference condition: 


RT CÉ, Poo ref 
Um = UM o + In | <E 2 16 
M M.o + ZF ( Bae (16) 


Um,o is the standard half cell potential of methanol oxidation, 
Pco,,ref is chosen 1 atm, Cy re is set 0.5 mol dm~°, that is, the 
maximum concentration used in the present study. It should be 
mentioned that the concentrations used in Eq. (16) should be in 
mol dm~>. The proton concentration Cy+ is available if the SPE 
conductivity Km is known according to the following relationship 


[33]: 


F2 
Km = pre Cat (17) 


(23) 


2.5. Current conservation in the anode catalyst layer 
The cell current density is calculated by integration of the 


local transfer current density across the thickness of anode cat- 
alyst layer. That is 


La 
Teen = f jdz’ (24) 
0 


Finally, 


| k 
keg = nM F kpett (co ( zat) 
a 
= h : L (CÌ + Ch) (25) 
cosec pea hs 7+Ch 
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2.6. Mass conservation of oxygen in the cathode catalyst 
layer 


Mass balance of oxygen is similar to that of methanol in the 


anode catalyst layer: 


2 
eff d Co 
c dz? 


—k-Co =0 (26) 


DË is the effective diffusivity of dissolved oxygen within the 
cathode catalyst layer, and ke is the rate constant dependent on 
the cathode over-potential: 


(27) 


: —ac Fe 
ke = AcloO, ref 


——— exp 

no FCo,ref RT 
Qo and igo set are the cathodic transfer coefficient and the refer- 
ence exchange current density, respectively. nc is defined 


Nc = Ve — Pe — Uo (28) 
where 

RT 
Uo = V0.0 + gF In( Po, ref CHy+ ref) (29) 


Co is to be determined with suitable boundary conditions. 
Let z” = z — (Lp + La + Lm + Le), the oxygen dissolution at 
CL/gas diffuser interface is described 


C 
Co= 75 atz" =0 (30) 


(0) 
where Cog is the oxygen gas concentration at the interface and 
Ko is the solubility of oxygen in the flooded cathode catalyst 
layer. At the membrane/cathode catalyst layer interface, there is 
no penetration of oxygen into the membrane: 


dCo 
g =O ee (31) 


Aqueous oxygen distribution within the cathode catalyst layer 
is determined 


1 + tanh ke Lo | tanh Ke 32 
+ tan pet Le | ta D? (32) 
C 


Partial current density contributed from the reduction of oxygen 
is 


x 


0 
Io = noF f keCo dz” 
ie 


Co | k 
aa / ff “V8 £ 


On the other hand, a mass conservation of oxygen in the cathode 
gas diffuser gives 


eff 


D 
d (Cog,b — Cog) (34) 


Io = noF 
o oF ia 


in which Cog,» is the oxygen gas concentration in the cathode 
flow channel, pf the effective gas diffusivity in the gas diffuser. 
Eq. (33) is further simplified to 


noF((y keD Cog b)/Ko) tanh(y/ke/ DFL.) 


(35) 
1+ (keD La/Ko DS) tanh(\/ke/ DF Le) 


o= 


2.7. Current conservation in the cathode catalyst layer 


The resultant current density in the cathode catalyst layer is 
the difference between the parasitic oxygen reduction current 
density and oxidation current density of methanol that crossed, 
or, the leaking current density. 


Teen = lo — leak (36) 


The cell voltage Vey is calculated by subtracting the activation 
losses in both anode and cathode (na and |7¢|), and the polymer 
electrolyte ohmic resistance (Lm {ecel /Km) from the theoretical 
thermodynamic cell voltage (Uo — Um): 

LL cen 


Veet! = Uo — Um + Ne — Na — ——— (37) 
Km 


2.8. Numerical procedure 


Solution procedure is based on the potentiostatic mode, na, Nc 
and Keen are calculated using Newton-Raphson method [34] sug- 
gested by Eqs. (25), (36) and (37) for a specified Vee. Numerical 
derivatives are used to calculate the jacobian matrix element Ak, 
at kth iteration [34]. That is, 


Fik + £64) — Fiat) 
Gi k k k 
(xj + ex) =x} 


be oe 
7 OX; | yk 


(38) 


F; represents the non-linear equation i to be solved, while x; 
represents Na, Nc OF Icen, 54; the Kronecker delta, and £ is a step 
size, 1076 used here. 


2.9. Parameter estimation 


The model is used to retrieve useful kinetic parameters by 
non-linear regression with available experimental data. In the 
present study, Da, Km and iom ref are estimated to correlate the 
experimental polarization curves [23] with the proposed model 
using Levenberg-Marquardt algorithm [35]. The minimization 
function is defined 


Fmin = > Uexp,e(Eexp,t) — Isim, (Eexp, i) (39) 
k 


in which the difference of the experimental and simulated cell 
current densities corresponding to a specified cell voltage is to be 
minimized in the sum of least squares. Fixed or explicit expres- 
sions of thermodynamic, kinetic, and transport parameters used 
are listed in Table 1 unless otherwise specified. For the case study 
employed, determination of the anode catalyst layer thickness 
La and the cathode catalyst layer thickness Le is given as follows. 
For the anode catalyst loading mpry at the 1:1 atomic ratio of 
PtRu alloy, the weight percentages of Pt and Ru in the alloy 
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Table 1 

Fixed parameters used in the model unless otherwise specified 

Parameter Expression References 

Aa 52 x 104 cm? g7! [36] 

Ac 52 x 10* cm? g! Chosen the same as A, 
Da 0.5035 x 1074 cm? s7! Estimated in the present study 
pet Dye! [33] 

Dp Da Chosen the same as Da 
pg" Dret” [33] 

D: 7.4 x 10-8(./@gMpT7/ V$ up) [38] 

Dg Deel” [33] 

Da 2.745 x 1074 / Pe(T/-/TeaTep) (Pea Peo)? (Tea To)? 1 /Ma) + 7M) [39] 

pF Daeg’ [33] 

Dy+ 1 x 1074 cm? s7! [40] 

Dm 4.9 x 1076 exp(2436((1/333) — (1/T))) [18] 

ioM.ref 0.7630 x 1077 A cm~? Estimated in the present study 
ip0 ref 0.3189 x 1077 A cm~? [37] 

kp 1.17 x 10!° exp(—19098/T) [20] 

Ky 1.25 [25] 

Ky 0.8 [23,25] 

Ko 1/RgT exp(—666/T + 14.1) [36] 

Lp 0.026 cm Chosen the same as La 
La 0.026 cm [33] 

Lm 0.023 cm [33] 

Ma 32g mol! [39] 

Mp 28g mol! [39] 

Mpg 18g mol! [39] 

Pea 49.7 atm [39] 

Po 33.5 atm [39] 

Tea 154.4 atm [39] 

Teh 126.2 atm [39] 

UMo 0.03 V [23] 

Uoo 1.23 V [37] 

Va 14.8 cm? mol! [39] 

Qa 1 Chosen the same as œe 
Qe 1 [36,37] 

Km 0.6842 x 107! scm! Estimated in the present study 
HB 1.45 cP [38] 

4 45 cm? mol! [31] 

Ec 0.3 [41] 

Ed 0.4 [37] 

Eb 0.4 [37] 

Ea 0.3 [41] 

Dg 2.26 [38] 


are calculated to be 0.659 and 0.341, respectively. For a 40 wt- 
% PtRu in carbon supported PtRu/C, the thickness of anode 
catalyst layer is determined 


= (0.659mpRu/ prt) + (0.341 mpRu/PRu) + (1 — 0.4)mpiru/0.40c) 


The same argument is applied to the cathode catalyst layer. For 
a 40 wt-% Pt/C with a loading mp;, Le can be determined 


a 
l — £a 


PPt, PRu and pc are the densities of Pt, Ru and carbon. The void 
fraction in the catalyst layer is denoted £a. The active area per 
unit volume a, is deduced if the active area per unit mass of 
PtRu, Ag, is available. 


A; 
pe aMPtRu (41) 
La 


(40) 


_ (mpr/pp) + (A = 0-Hmp/0.4pc) 


L 
e l — £c 


(42) 


€c is the void fraction in the cathode catalyst layer. The active 
area per unit volume a, is expressed 


Acmpt 
Le 


(43) 


de = 


K.-M. Yin / Journal of Power Sources 167 (2007) 420-429 427 
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0.00 0.02 004 006 008 010 O12 0.14 0.16 
| /Acm? 


cell 


Fig. 2. Comparison of theoretical (lines) and experimental (symbols) polariza- 
tion curves at different fed methanol concentrations. Experimental data is from 
Garcia et al. [23]. Anode catalyst: 40% PtRu/C, atomic ratio Pt:Ru= 1:1, PtRu 
loading: 3 mg cm; cathode catalyst: 40% Pt/C, Pt loading: 1 mg cm~? (E-TEK 
Inc.). Membrane: Nafion® 117. The cell and feed flow temperatures were set 
70°C. Estimated parameters from the model: io ret = 0.7630 x 1078 Acm~?, 
Da =0.5035 x 10~4 cm? s~!, km = 0.6842 x 107! S cm™!. 


3. Results and discussion 


The experimental polarization data in Ref. [23] and model 
predictions in the present study at various fed concentrations of 
methanol are depicted in Fig. 2. As illustrated in the plot, the 
limiting behavior occurs much earlier when lower methanol con- 
centration is used, that is, it initiates at larger cell voltages. The 
theoretical limiting current densities can be calculated simply by 
setting methanol concentration zero at the anode backing/anode 
catalyst layer interface, so that 


lim = nM ppe & (44) 
Ly 
The limiting current densities calculated from Eq. (44) are con- 
sistent with the numerically predicted values in Fig. 2 at limiting 
discharge rates. Note the effects possible micro-convection and 
gas void within the diffusion layer on the cell discharge, induced 
by the CO2 generation, are included implicitly in the estimated 
diffusivity of methanol in the backing layer Dp (=D, in Table 1), 
as discussed in Assumption 2 in Section 2. The predicted dis- 
charge curves near the cell open circuit voltage (OCV) exhibits 
a higher cell voltage when the methanol concentration is lower. 
This is a consequence of the reduced cathode mixed potential 
due to the crossover of methanol. Significant voltage loss occurs 
at the cathode as more methanol crossed due to a higher anode 
feed concentration. This point will be discussed in more details 
latter. Fig. 3 shows the decrease of crossover current density with 
the cell current density. As expected, maximum eak appears at 
the open circuit voltage in which mass transfer rate of methanol 
is the highest. As the cell discharge rate approaches the limit- 
ing value, where nearly no concentration of methanol available 
in the anode catalyst layer, the crossover is diminished accord- 
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Fig. 3. Leaking (crossover) current density vs. cell current density at various 
methanol concentrations. Key is the same as in Fig. 2. 


ingly. The corresponding methanol concentrations within the 
anode catalyst layer and membrane, and oxygen concentration 
within the cathode catalyst layer are profiled in Fig. 4 at the cell 
potential 0.5 V. Non-linear concentration profiles appear within 
the catalyst layers. The discontinuity of methanol concentration 
at the anode catalyst layer/membrane interface is governed by 
the phase equilibrium constraint [23,25]. It is shown at the low 
methanol concentration of 0.05 M limiting behavior emerges as 
earlier as 0.5 V (also referred to Fig. 2). Oxygen concentration 
declines rapidly in the cathode catalyst layer. Deepest oxygen 
penetration path appears at the methanol concentration 0.05 M 
since lower cell discharge rate allows more residence time for 
the diffusion of oxidant through the medium. 
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Fig. 4. Simulated methanol concentration profiles within the anode catalyst 
layer, membrane, and oxygen concentration profiles in the cathode catalyst layer 
at Vece =0.5 V. Key is the same as in Fig. 2. 
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Fig. 5. Anode and cathode potential responses as a function of cell current 
density at [MeOH] =0.05 and 0.5 M. 


It is informative to inspect the respective anodic and cathodic 
potential responses as a function of cell current density, as indi- 
cated in Fig. 5. The open circuit potentials of methanol oxidation 
at anode Uy, and oxygen reduction at cathode Uo are plotted as 
the references. The cathodic potential, Ve — c, and the anodic 
potential, Va — ®a, are displayed respectively. Note the voltage 
difference between these two is equal to the cell voltage Veen 
plus the ohmic loss Lin Kecel /Km according to Eq. (45), as derived 
by the combination of Eqs. (15), (28) and (37): 


ml, cell 


= (Ve — Bc) — (Va — Pa) (45) 


L 
Veell + 


m 


It is seen the anodic polarization is higher when lower methanol 
feed concentration (0.05 M) is used, and the anode mass transfer 
limitation is arrived much earlier than that at high feed concen- 
tration (0.5 M). As for the potential characteristics at cathode, 
higher voltage gain is observed at low than that at high methanol 
feed as most prominent near the open circuit voltage. This is the 
consequence of two factors. Firstly, methanol crossover is most 
severe at low discharge rate and at high anode concentrations 
as illustrated in Fig. 3. Methanol oxidization at cathode con- 
sumes a portion of the cathode oxygen reduction current so as 
inducing a less available cathode potential that should be oth- 
erwise higher. Second, the fact that cell current density limited 
by the low methanol concentration implies a less requirement 
of cathode reduction rate. A lower cathode over-potential |7¢| is 
expected as a result. It is also shown in Fig. 5 that the available 
cathode potential is generally large and sufficient for a com- 
plete oxidation of crossed methanol. That is, an instantaneous 
consumption of the permeated methanol is anticipated. 

Fig. 6 depicts the influence of cathode oxygen feed concen- 
tration. It is obviously that higher cell OCV and an improved cell 
discharge behavior are predicted as pure oxygen is used. Mass 
transfer limitation is still governed by the anode electrochemical 
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Fig. 6. Effect of oxygen feed composition on the polarization curve. 


reaction whether air or oxygen is fed. The corresponding oxygen 
concentration profiles in the cathode catalyst layer (not shown) 
indicate more oxygen is preserved in this layer when oxygen 
instead of air is fed. On the other hand, methanol concentration 
depletion is more significant in the anode catalyst layer when 
oxygen is used. 

Fig. 7 illustrates the /-V characteristics as the cathode pres- 
sure is increased. Higher pressure in the cathode suggests a 
prohibition of methanol crossover so that a less mixed poten- 
tial effect is expected. Further, higher oxygen concentration is 
available for the reduction reaction when cathode is pressur- 
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Fig. 7. Cathode air pressure effect on the polarization curve of DMFC. 
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Fig. 8. Concentration profiles of methanol and oxygen within the MEA at 
Veet! =0.5 V, Pe — Pa =0, 1, and 2 atm. 


ized. Fig. 8 depicts the detailed concentration distributions at 
Veen = 0.5 V. As expected, higher oxygen concentration is main- 
tained within the cathode catalyst layer with raised cathode 
pressure. On the other hand, more methanol consumed as more 
oxygen is available in the cathode. The inhibition of methanol 
crossover is clear as a significant depressed concentration profile 
is observed in the membrane phase when the cathode pressure 
is increased. 


4. Conclusion 


An algebraic one-dimensional mathematical model on the 
direct methanol fuel cell is developed. Mass transport of 
methanol or oxygen through the various media of the membrane- 
electrode-assembly can be correlated well with the associated 
non-linear electrochemical kinetics within the catalyst layers. 
Methanol crossover effect is accounted by the mixed potential 
theory and justified by the predicted anodic and cathodic poten- 
tial responses. With the incorporation of a non-linear parameter 
estimation scheme, this model can be used for a better design 
on the MEA fabrication of DMFC. 
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